
function [Q,varargout] = alphaVariance(yt,B,R,q,Q_0,T_0)


Y = yt(:,q+1:end);
X =[];
for qx=q:-1:1    
    X = [X;yt(:,qx:end-q+qx-1)];
end

nu = R\(Y - B*X);
[Q,sse] = drawVariance(nu,Q_0,T_0);

varargout(1) = {nu};
varargout(2) = {sse};


